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Model for the magnetoresistance and Hall coefficient of inhomogeneous graphene 



00 

o 
o 

O 

Q 

00 



Rakesh P. Tiwari and D. Stroud 

Department of Physics, 
Ohio State University, Columbus, OH 43210 

(Dated: December 9, 2008) 

We show that when bulk graphene breaks into n-type and p-type puddles, the in-plane resistivity 
becomes strongly field dependent in the presence of a perpendicular magnetic field, even if homoge- 
neous graphene has a field-independent resistivity. We calculate the longitudinal resistivity p xx and 
Hall resistivity p xy as a function of field for this system, using the effective-medium approximation. 
The conductivity tensors of the individual puddles are calculated using a Boltzmann approach suit- 
able for the band structure of graphene near the Dirac points. The resulting resistivity agrees well 
with experiment, provided that the relaxation time is weakly field-dependent. The calculated Hall 
resistivity has the sign of the majority carrier and vanishes when there are equal number of n and 
p type puddles. 
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Graphene is a two-dimensional form of carbon with 
a hexagonal crystal structure like that of a single layer 
of graphite. Because of this structure, it has the band 
structure of a semimetal: the Fermi energy Ep of neu- 
tral graphene lies at a "Dirac point," where the electronic 
density of states n(Ep) — 0. There are two inequivalcnt 
Dirac points located at different Bloch vectors ko and ki . 
Near the Dirac points, the bands are linear functions of 
the components of k — k and k— ki, and n(E) is propor- 
tional to \E — Ep\. Because of this unusual band struc- 
ture, the quasiparticle Hamiltonian near the Dirac points 
is formally identical to that of massless Dirac fermions, a 
feature which is responsible for part of the recent interest 
in graphene. 

Graphene also has striking transport properties. For 
example, experiments have observed finite conductivity 
for all values of Ep, whether above or below the Dirac: 
point[l[, with a minimum conductivity typically ~ 4e 2 jh. 
However, some workers have suggested that this mini- 
mum could have much smaller or larger fH values than 
4e 2 //i. It has been proposed that the existence of a finite 
conductivity even at the charge neutrality point might 
be a result of local potential fluctuations, which could 
cause a homogeneous neutral graphene sheet to break 
up into "puddles" of electron-rich (n-type) and hole-rich 
(p-type) character [U. These puddles have, in fact, been 
unambiguously observed in experiments using scanning 
tunneling microscopy Q. 

Recently measurements of the magnetic-ficld- 
dependent longitudinal and Hall resistivity p xx and 
p xy measurements have been reported 0- Pxx was found 
to increase by nearly tenfold with increasing magnetic 
field perpendicular to the graphene film, followed by 
an apparent saturation at sufficiently strong magnetic 
field. These authors found that the magnetoresistance 
was inconsistent with a two-fluid model of transport 
by n-type and p-type charge carriers in a homogeneous 
sheet of graphene, and suggested that it might agree 
better with a model describing the film as a mixture of 



n-type and p-type puddles. However, they were able to 
obtain close agreement between experiment and theory 
only by assuming an ad hoc empirical form for the 
magnetoresistance. 

In this Rapid Communication, we present a simple 
model for the magnetoresistance and Hall coefficient of 
graphene, based on the effective-medium approximation 
(EMA) in a transverse magnetic field. Such a model 
is reasonable if the n-type and p-type puddles are dis- 
tributed randomly, as appears to be the case in Ref. 0]. 
Our results show that when the area fractions of n-type 
and p-type puddles are exactly equal, p xx varies exactly 
linearly with field. At other puddle fractions, it satu- 
rates, in agreement with experiment. We find that we 
can obtain excellent agreement with the observed behav- 
ior of p X x(B) if we assume an n-type area fraction /„ 
satisfying |/„ — 1/2 1 ~ 0.07 and reasonable values for av- 
erage carrier density and transport relaxation time. We 
also make predictions about the Hall resistivity p xy (B). 

We consider magnetotransport in a single layer of 
graphene subject to a magnetic field B = Bz perpen- 
dicular to the graphene layer. We assume that, because 
of a random potential due to charges in the substrate 
or some other cause, the graphene layer has broken up 
into a mixture of n-type and p-type puddles, having area 
fractions /„ and f p = 1 — /„. We also assume that each 
of the puddles is large enough to be described by its own 
magnetoconductivity tensor a n or a p . In practice, this 
assumption means that the puddle dimensions are larger 
than a typical carrier mean free path. This condition 
may not always be satisfied in practice, in which case 
the results below might need to be modified. 

We assume that a n and a p are both given by the usual 
free-electron (or free-hole) forms, suitably modified to ac- 
count for the linear dispersion relations of the electrons 
and holes near the Dirac point. Thus, for the 2x2 con- 
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ductivity tensor in the xy plane we write 
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We denote the charge carrier densities in the n-type and 
p-type puddles by n and p, and the corresponding re- 
laxation times by r n and t p . We also assume that all 
the puddles have the same density of charge carriers, so 
that n = p, and that the relaxation times t„ = r p = r. 
With these assumptions, the zero-field conductivities ao 
of the n-type and p-type puddles are equal. We can also 
define zero-field mobilities p n and fi p by er n = ne/i n and 
er p = pe/ip, where e is the magnitude of the electronic 
charge; with the above assumptions, these mobilities are 
also equal. 

Both the zero-field conductivity ao of the puddles and 
the cyclotron frequency uj c are modified from their usual 
free-electron values because of the linear dispersion rela- 
tions near the Dirac point. The result for ao at temper- 
ature T = (fully degenerate limit) is 



f n 7^ 1/2, corresponding to a net doping. This would 
correspond to a graphene sheet biased by a suitable gate 
voltage. 

A reasonable way of calculating a e for tensor conduc- 
tivities is provided by the effective-medium approxima- 
tion (EMA) 9]. In this approach, the electric fields and 
currents within the inhomogeneous graphene sheet are 
calculated as if the n-type and p-type puddles are com- 
pact and approximately circular, and are embedded in 
an effective medium whose conductivity is calculated self 
consistently d 03- The EMA is, in fact, exact at f = 
1/2, provided fi p = [i„ [lof. For tensor conductivities, 
the defining equation for the EMA is 
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Here bai = Oi — a e , I is the 2x2 unit matrix, and, for 
the planar geometry considered, T = —I/(2a etXX ) is the 
depolarization tensor. This matrix equation reduces to 
two coupled scalar algebraic equations for the two inde- 
pendent components of a e {a e>xx and a e ^ xy ), which are 
easily solved numerically. The other two components are 
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where vp is the Fermi velocity. This form is obtained 
from the usual solution of the Boltzmann equation for a 
degenerate Fermi gas@, which gives 

a xx = [2e 2 T/(27r 2 )] / Sk'tT 2 {dE{]s!)/dk' x f 5 (E(k') - E F ) 



. (4) 

for the conductivity. Here, k' is the two-dimensional 
wave vector, measured relative to one of the Dirac points, 
and E(k') = hvp |k'| is the energy relative to the Dirac 
point. We use kp = y/im, which takes into account 
the two valleys near the two inequivalent Dirac points 
in the graphene band structure, and we have included 
an extra factor of 2 in eq. for the same reason. 
The cyclotron frequency lo c is readily obtained from the 
semiclassical equation of motion fik = evk x B, where 
Vk' = H~ 1 'V\ c /E(]s.') ) as applied to a band with the disper- 
sion relation E(k') — vpH\]c'\; the result is (in SI units), 
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Next, we calculate the effective conductivity tensor a e 
of a graphene sheet which has broken up into n-type and 
p-type puddles. If /„ = 1/2, this would correspond to 
the case where the net charge carrier density is zero, cor- 
responding to a neutral graphene sheet which would, if 
homogeneous, have its Fermi energy at the Dirac point. 
However, it is also possible to have a graphene sheet with 



resistivity tensor is then obtained by inverting the ma- 
trix cr e , so that p &iXX = p et yy = a e xx I {a^ xx + o\. xy ) and 

Pe,xy = ~Pe,yx = ~ cr e,xy / { cr e,xx ^e.xj/) - 

In order to compare this model to experiment [7[ , we 
need the values of vp, n, r (or equivalently, fi), and 
f n . From the band structure of graphene, vp ~ 10 6 
m /secQ. In fact, a value of 1.1 x 10 6 m/sec has been 
inferred from measurements of the Landau level splitting 
in graphene [1], and we use this value in the calculations 
below. Also, the measured value of the zero-field resis- 
tivity is CTq -1 = po — 0.125/i/e 2 . Given this value, eq. (3) 
provides one condition satisfied by the two parameters n 
and t. We then choose n, r, and f n so as to best fit the 
measured p X x{B) at B = 8T, and to yield lo c t = 3. IB, 
where B is the magnetic field in T, as reported in Ref. 0]. 
This procedure gives n~6x 10 14 m~ 2 and /„ ~ 0.43. 
We find that the best agreement with the resistivity is 
given at high fields by lu c t ~ 2.3 B, and at low fields by 
lo c t ~ 3.1 B, indicating a weakly field-dependent r. The 
value of n is close to measured value quoted in Ref. Q- 

The calculated results for p X x(B) are shown in Fig. 
1, using these parameters. As can be seen, the fit to 
the experimental data is excellent over most of the field 
range, using uj c t ~ 2.3B, and at low fields using uj c t ~ 
3. IB. The fit, especially at high fields, is also superior to 
the two-fluid model discussed (and found inadequate) in 
Ref. 0. The fit to this puddle model would be nearly 
perfect over the entire range of B studied experimentally, 
if t varied by about 30% as a function of B. The results 
for p xx are independent of the sign of the charge and are 
thus unchanged if /„ — * 1 — /„. 

In Fig. 2, we show the corresponding results for p xy (B). 
In this case, we use a field-independent r (corresponding 
to lj c t = 2.3 B). We show results for /„ = 0.43, 0.57, and 
0.5. p xy (B) for f n = 0.43 is equal and opposite to that 
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for /„ = 0.57. In both cases, p X y varies roughly linearly 
with B for B greater than about IT. At f n = 0.5, p xy = 
for all B. Within the present model, this lattter result is 
exact, and not restricted to the EMAflpj. 

Fig. 3 shows p xx (B,fn) versus f n for several values of 
B. We use the EM A and the same parameters as in Figs. 
1 and 2 (with lu c t — 2.3 B). p xx saturates for all val- 
ues of f n except /„ = 1/2, for which it increases linearly 
with B. Once again, this linearity is exact, and not re- 
stricted to the EM A [To1| . In Fig. 4, we show p xy (B,f n ) 
versus f n for several values of B. As can be seen, p xy 
changes sign at /„ = 1/2, and approaches a constant as 
f n approaches either 1 or 0. The magnitude of the slope 
[dp X y/df n ]f n= i/2 increases with increasing B, so that, at 
large B, the Hall resistivity is very close to that of the 
majority charge carrier. 

The present model agrees well with the measured val- 
ues of p xx {B) in graphene. However, it is based on cer- 
tain assumptions whose validity for graphene we now dis- 
cuss. One assumption is that graphene can be treated as 
a macroscopically inhomogeneous assembly of puddles, 
each with its own conductivity tensor. The scanning tun- 
neling microscope images shown in Ref. 0] suggest that 
the carrier density varies appreciably over a distance of 
perhaps 0.2/i. Using the above density estimates, the 
number of charge carriers in a puddle of linear dimen- 
sion 0.2/Lt would be ~ 25. This size is rather small to be 
treated macroscopically. On the other hand, a more rea- 
sonable definition of a "puddle" might be a region where 
the charge carriers were all of one sign. Judging from 
the images, a typical linear dimension of such a region 
would be larger than 0.2/i - perhaps 0.5 - lp, and would 
contain ~ 500 carriers. This is probably large enough to 
describe each puddle by its own macroscopic conductiv- 
ity, provided that the mean free path A is less than lp. 
Even if A were larger than this, one could still use the 
semiclassical model with a size-limited A of magnitude 
A ~ a, where a is the linear dimension of the puddle. 
This would give t ~ a/vF ~ 4.5 x 10 -13 sec. In short, 
treatment of the puddle mixture macroscopically is prob- 
ably appropriate in the case of some disordered samples 
of graphene, and the results seem to agree with experi- 
ment. 

Another point, as can be seen from Fig. 3, is that the 
Pxx{B, f n ) saturates (approaches a finite limit at large B) 
only if /„ 7^ 1/2, i. e., if there is a net charge imbalance 
produced by a suitable gate voltage. Such saturation 
seems to be observed in experiments^. But, as can be 
seen from Figs. 2 and 4, a charge imbalance would lead 
to a nonzero p xy (B). It would be of interest if p xy (B) 
could be measured and compared to the values needed 
for the present model to agree with experiment. 

Thirdly, the present model treats the electron dynam- 
ics semiclassically, and thus does not take into account 
the quantum Hall effect (QHE), which is seen at suf- 
ficiently high fields (1 03]. Typically, the QHE will be- 
come visible when the spacing between the Landau levels 
is large compared to fc^T. This can occur even at room 



temperature in graphene [12|. If the QHE becomes im- 
portant, the present semiclassical model would need to 
be modified. 

Finally, we obtain the best fits to experiment if we as- 
sume a weakly magnetic-field dependent relaxation time, 
as described above. Such field-dependence could be rea- 
sonable, but it would be useful to have a model which 
explicitly produces a magnetic field-dependent r. 

We find that our results are quite insensitive to slight 
changes in the parameters or other features of the model. 
For example Fig. 1 suggests that p xx changes only 
slightly, but not dramatically, when t is varied by ~ 30%. 
Also, we have recalculated p etXX (B,T) without the as- 
sumption that the electrons and holes have equal mo- 
bilities. Even if the mobilities are different, we find 
that p xx {B) still varies linearly with B at /„ = 0.5 and 
saturates at other values of /. Another change in our 
model is suggested by that fact that the carrier density 
in graphene must be a continuous function of position, 
rather than being simply bimodal as postulated in our 
model. To check the effects of a non-bimodal distribu- 
tion, we have repeated our calculations assuming four 
types of puddles, two n-type and two p-type, with two dif- 
ferent densities each of electrons and holes. Once again, 
we find that the resulting p xx (B,f n ) depends primarily 
on /„ and r, and not on the presence of two types of n 
and of p puddles. Finally, we have considered the case 
of a three-component composite, made up of n-type, p- 
type, and insulating regions. Here, once again, we find 
that, for a small insulating areal faction (~ 0.1), we ob- 
tain linear magnetoresistance if /„ = f v and saturating 
magnetoresistance otherwise, similar to the case of no in- 
sulating regions. Thus, our results are not much affected 
by small modifications in our model. 

To summarize, we have calculated the magnetoresis- 
tance and Hall resistivity for a semiclassical model of 
graphene, on the assumption that it is a mixture of n- 
type and p-type puddles, and using the correct form of 
the band structure near the Dirac points. The resulting 
magnetoresistance is in good agreement with experiment, 
provided that the areal fractions of n and p-type pud- 
dles are slightly different and that the relaxation time is 
weakly magnetic-field dependent. Further confirmation 
of the model could be obtained if the measured Hall resis- 
tivity were compared to that computed from this model. 
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FIG. 2: (Color online.) Calculated p xy (B, /„), as obtained in 
the EMA, using fj, = 2.3T -1 and /„ = 0.431, 0.569, and 0.5. 



(a) 



— u. = 3.1T 1 ,f =0.4287 or0.5713 

~ n 

-- ii = 2.3T 1 ,f =0.431 or 0.569 

~ n 

o Experimental data points 





B = 2T 


2 — 


B=4T 


3 — 


B = 6T 


4 — 


B = 8T 


5 — 


B = 10 T 



( b) 




FIG. 1: (Color online.) p X x(B,f n ) as a function of B (in 
T), with two different assumptions about the mobility, In 
both cases, we assume that the electron and hole mobili- 
ties are equal. Solid (red) line: calculated results, assum- 
ing /i = lu c t/B = 3.1T" 1 and /„ = 0.429 (or 0.571).. 
Dashed (black) line: calculated results with p = 2.3T -1 and 
/„ = 0.431 (or 0.569). Open circles are experimental data 
from Ref . [7j . Lower panel is a blowup of the calculations and 
data from the upper panel. 
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